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The cross section for muon pair productions by electrons scattering 
O \ over photons, <Jmpp, is calculated analytically in the leading order. It 
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is pointed out that for the center-of-mass energy range, s > 5m^, the 

cross section, aupp is less than 1 /ib. The differential energy spectrum 

^ ■ for either of the resulting muons is given for the purpose of high-energy 

, neutrino astronomy. An implication of our result for a recent sugges- 

00 '. tion concerning the high-energy cosmic neutrino generation through this 

_ , ' muon pair is discussed. 
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00 

T— I ' High-energy neutrino astronomy is now a rapidly developing field This entails 

^ I the identification of possible sources for high-energy cosmic neutrinos. In particular, non 
^ I hadronic processes taking place in cosmos remain an interesting possible source for high- 
^ ! energy neutrinos. 

Ph! In this context, Kusenko and Postma have recently suggested that, for a class of topo- 

Ph. logical defects that may produce ultrahigh-energy photons at high red shift, the dominance 
^ ■ of muon pair production (MPP) in e~7 — > e~fi~^fi^ over triplet pair production (TPP) 
' in e~7 e~e^e~ for s > 5m^ enables the MPP process to be an efficient mechanism 
^ ■ for generating high-energy cosmic neutrinos 0. The electrons in the initial state of the 
' above processes are considered as originating from the electromagnetic cascade generated 
by the ultrahigh-energy photons scattering over the cosmic microwave background photons 
present at high red shift. Subsequently, assuming that this ultrahigh-energy photon source 
at high red shift is responsible for the observed flux of ultrahigh-energy cosmic rays, it 
was pointed out that the decays of muons in MPP give rise to high-energy cosmic neu- 
trino flux -E^^^ typically peaking at E^, ~ 10^^ eV 0. Concerning the cross sections of the 
above processes, we note that, although TPP process has been thoroughly studied before 
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1^, there is still no explicit estimate for (Tmpp for s > 5m^ to our knowledge. In order 
to verify whether MPP indeed dominates over TPP for s > 5m^, we analytically calculate 
the (Jmpp for < s < 20 m^. In this Rapid Communication, we present some details 

of the calculation and compare our result with the value quoted in Refs. Our con- 

clusion is that, within the above s range, the cross section aMPP we have obtained is at 
least three orders of magnitude smaller than the one quoted and used in Refs. 0,^. This 
implies that MPP can not be the dominating high-energy neutrino generating process as 
suggested in Refs. 0J§]. This conclusion is based upon the value of the ratio R defined as 
R ~ aMPP / iVTPpO'TPp) [i], where r]Tpp{s) ~ 1.768(s/mg)~°-''^ is the inelasticity for TPP 
and (Ttpp — a^m~^ [3.11 In (s/rrig) — 8.07], both for s ^ 0]. The rjTpp is basically 
the average fraction of the incident energy carried by the final state positron. The original 
estimate of Refs. gives R ~ 10^, which favors MPP process as the dominating cosmic 
high-energy neutrino generating process. Namely the electron energy attenuation length due 
to TPP process is much longer than the interaction length of MPP process. However, since 
the correct value for (Tmpp within the energy range 5m^ < s < 20 is three orders of 
magnitude smaller, the ratio R becomes less than 1. Therefore MPP is no longer an effective 
mechanism for generating high-energy cosmic neutrinos at the high red shift. 

There are 8 Feynman diagrams contributing to the MPP process in the leading order. 
These are shown in Fig. 0. Among these 8 diagrams, 4 diagrams contain a Z-boson exchange, 
which can be disregarded for the range of s under discussion. Among the remaining 4 
diagrams, one can also disregard diagrams (c) and (d) due to the inflow of large invariant 
energies into the electron [diagram (c)] or photon [diagram (d)] propagators. We therefore 
only perform analytic calculations for diagrams (a) and (b) in Fig. while evaluate the rest 
of the diagrams numerically. We have applied the equivalent photon approximation to 
compute diagrams (a) and (b). We write the MPP cross section in the following convolution 

ctmpp = [ dxf^/eix)ayy^f,+^-{s = xs), (1) 

J4m2 /s 

where f-y/e{x) = (Q;/27r)[(H- (1 — x)^)/a;] ln(s/mg) is the probability of finding a photon from 
the incident electron with an energy fraction x = E^/E^. The differential cross section as a 
function of the outgoing muon energy i?^, in head on collisions, is given by 



dCTA-J 



MPP 
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+ tanh ^(f) 



(2) 



where y = E^/m^ with ymm = 1 and ymax = y/s/^nifj,; v is the velocity of the outgoing muon, 
which is related to ?/ by f = yT^^^. Apart from the trivial prefactor, the first line in Eq. 
arises from the distribution function f^/e{x) while the second line describes the cross section 
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for 77 jj,'^ ii~ . The total cross section aMPP, obtained by performing the y integration, 
depends only on the center-of-mass energy y/s. For s = 5m^, aMPP — 4 ■ 10~^/ib, while 
o'MPP increases to about 0.1 yub for s = 20 m^. These cross-section values are at least 3 
orders of magnitudes smaller than those of Refs. 0,0], which give (Tmpp between 0.1 and 
1 mb in the above energy range. The last equation gives the differential energy spectrum 
of either of the muons produced in MPP. In the large s limit, the (Tmpp given by Eq. 
behaves as 

2q/3 ( s \ 

.MPpW.^ln(2)lny . (3) 

The inelasticity for MPP (in the center-of-mass frame) is defined as 

VMPP{s) = I dE, { ^\ (4) 



Using Eqs. and (|^), we find that the average fraction of incident energy carried by either 
of the muons in MPP, in the large s limit, behaves as 

/ s -0.5 

nMPp{s) ~ 3.44 y-^j , (5) 

whereas it is close to 0.48 near the threshold for the MPP process and approaches to 0.34 
for s ~ 20 m^, . 

We have used a package for computations in high energy physics, CompHEP (version 
33.22), to check our results for ctmpp [H- The comparison for 5m^ < s < 20 is shown in 
Fig. ^. We note that the plot generated by CompHEP is a result of computing all 8 diagrams. 
From Eq. (H), we also note that, for s ^ 5m^, the MPP cross section aupp < 1 /ib. We 
have confirmed this remark using CompHEP as well. We have made a further check on our 
O'MPP by generating the amphtude of diagrams (a) and (b) symbolically using the package 
FeynArts 3 0. In this procedure, we do not use the equivalent photon approximation. We 
then obtained <Jmpp by numerically performing the phase space integration. We have found 
a rather good agreement (within few percent) between the (Tmpp obtained in this way and 
that given by the equivalent photon approximation. Because of the rather lengthy algebraic 
expressions occurring in the above procedure, we are omitting further details of this check. 

To obtain a consistent estimate of R, it is desirable to compute otpp and rjTppm addition 
to O'MPP- This can be easily done in the large s limit. By replacing /i^ with in the final 
state, the same 8 diagrams shown in Fig. |I| also contribute to the TPP process in the 
leading order. The dominating diagrams are again the first two. In the equivalent photon 
approximation with s ^ m^, they give 

arppis) ^ 4 1^(2) In i^] , (6) 

mi \ mi 
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whereas 



Vtpp{s) ~ 3.44 




The result for (Ttpp can be easily inferred from Eq. by replacing there with ml and 
multiplying a symmetry factor 1/2 which takes into account the identical-particle effect in 
the TPP process. The above results for o"tpp and rjTpp agree well (within few percent) 
with the results quoted in We further evaluated riTpp{s) in the Lab frame also, which 
in the large s limit, for the two dominating diagrams, is ~ 0.25(s/mg)~°'^. In practice, 
one should use the Lab- frame t]tpp in the the definition of R. Numerically, the Lab-frame 
r]Tpp is approximately an order of magnitude smaller than its center-of- frame counterpart. 
This suppression is related to the kinematical factors present in the Lorentz boost from the 
center-of- mass frame to the Lab-frame. Using the Lab-frame rjTpp and the values for o"tpp 
and (Tmpp, we still obtain i? < 1 for 5m^ < s < 20 m^. 

In summary, the authors in Refs. use the value 0.1-1 mb for the cross section 
o'mpp in the energy range 5m^ < s < 20 m^. As a result, they have deduced that R ~ 
o'MPp/iVTPpO'TPp) ^ 1- However, as we have shown, the correct value for aupp obtained 
from Eq. yields R < 1. In particular, near the threshold for the MPP process, i.e., for 
s ~ 5m^, we have i? <^ 1. Therefore, based on this observation, we conclude that MPP 
can not be a dominating high-energy cosmic neutrino generating process as suggested in 
Refs. Furthermore, we have also derived analytic expressions (in the leading order) 

for differential energy spectrum and the inelasticity of the MPP process in the large s limit, 
which might be of relevance in some other contexts of high-energy neutrino astronomy. 

For completeness, let us add a final remark concerning the possibility of high-energy 
cosmic neutrino generation in an electromagnetic cascade at high red shift, z, through 77 
collisions. The process 77 can, in principle, generate high-energy cosmic neutrinos 

near the threshold for this process, namely, for ^ 2m^. The neutrinos are generated 
through the subsequent decays of the final-state muons for a rather small range of z values, 
typically, for z ~ 5 — 10. In this z range, the interaction length for the process 77 is 
smaller than the energy attenuation length dictated by the process 77 e~^e~ . Furthermore, 
this interaction length is also smaller than the horizon length, cH{z)^^, where H{z) is the 
Hubble constant in this z range. 
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diagram (a) 



diagram (b) 
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diagram (c) 



diagram (d) 
► 1 — 



diagram (g) 



diagram (e) 



diagram (h) 



diagram (f) 



FIG. 1. The Feynman diagrams contributing to MPP in the leading order. 
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FIG. 2. (TAfpp as a function of square of center-of-mass energy s (in units of m^). Sohd curve is 
obtained using Eq. Dashed curve is obtained using CompHEP (it includes contribution from 
all 8 diagrams, see [6]). 
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